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ABSTRACT 

Many lens surveys have hitherto used observations of large samples of background 
sources to select the small minority which are multiply imaged by lensing galaxies 
along the line of sight. Recently surveys such as SLACS and OLS have improved the 
efficiency of surveys by pre-selecting double-redshift systems from SDSS. We explore 
other ways to improve survey efficiency by optimum use of astrometric and morpholog- 
ical information in existing large-scale optical and radio surveys. The method exploits 
the small position differences between FIRST radio positions of lensed images and the 
SDSS lens galaxy positions, together with the marginal resolution of some larger grav- 
itational lens systems by the FIRST beam. We present results of a small pilot study 
with the VLA and MERLIN, and discuss the desirable criteria for future surveys. 

Key words: Gravitational lensing ~ surveys 



1 INTRODUCTION 



Over the 30 years since the discovery of the first gravita- 
tional lens system (Walsh et al. 1979) about 100 such sys- 
tems have been discovered. The deflection of light in these 
systems can be used to probe the gravitational potential 
of the lensing galaxy. They can therefore be used to inves- 
tigate mass distributions in galaxies in the 0.5-20kpc re- 
gion where the images form. This is the interesting region 
where the transition between the dark matter-dominated 
halo and baryon-dominated core occurs. Combined with op- 
tical spectroscopy, this is the only unambiguous measure- 
ment of galactic mass available at cosmologically interesting 
distances. There have been numerous studies of mass distri- 
butions in galaxies, particularly in cases where the system 
gives enough constraints for a good lens model or when op- 
tical spectroscopy is also available or both (e.g. Kochanek 
et al. 2001; Saha & Williams 2001; Cohn et al. 2001; Munoz 
et al. 2003; Treu & Koopmans 2004; Dye & Warren 2005, 
Koopmans et al. 2006). A further use of galaxy lensing is 
the calculation of the Hubble constant using the time de- 
lay method of Refsdal (1964); the Ho value and error bud- 
get of this determination is somewhat controversial (e.g. 
Schechter et al. 1997; Impey et al. 1998; Saha & Williams 
2001; Kochanek 2002; Koopmans et al. 2003) but it never- 
theless represents an important, intrinsically clean method 
of Ho determination on cosmological scales. 

Most early lens searches relied on systematic surveys 
of background objects in order to pick out the few cases 



with galaxies close enough to the line of sight to produce 
multiple imaging. In the radio, the MIT-Greenbank survey 
(Bennett et al. 1986; Lawrence et al. 1986; Hewitt et al. 
1988) found five systems, a survey for lensed radio lobes by 
Lehar et al. (2001) and Haarsma et al. (2005) found a few 
further lens systems and the JVAS/CLASS surveys (Pat- 
naik et al. 1992; King et al. 1999; Myers et al. 2003; Browne 
et al. 2003) found 22 lens systems. In a subsequent sur- 
vey using a similar approach to CLASS (Winn et al. 2001; 
Winn et al. 2002) found at least three lens systems at South- 
ern declinations. Optical searches have included systematic 
ground-based studies of optically-selected quasars (Cramp- 
ton, McClure & Fletcher 1992; Surdej et al. 1993; Jaunsen et 
al. 1995; Kochanek, Falco & Schild 1995) and similar stud- 
ies with the HST (BahcaU et al. 1992; Maoz et al. 1992). 
More recently, the Sloan Digital Sky Survey has been used 
for further lens searches (e.g. Inada et al. 2003). Oguri et 
al. (2006) describe an algorithm based solely on the SDSS 
using morphological and colour selection, developing earlier 
work by Pindor et al. (2003). Pindor et al. use morphologi- 
cal selection by comparing the of single-component and 
multi-component PSF fits to SDSS images, and claim to be 
able to detect lens systems with separations as small as 0'.'7 
and with flux ratios as much as 10:1. 

Most searches are inefficient because typically many 
hundreds of background sources must be searched in order 
to find one gravitational lens, and it can be difficult to find 
an efficient filter which will select all the lenses while re- 
jecting false positives. The exception is the Haarsma et al. 
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(2005) survey which is discussed below, and, increasingly, 
the SDSS morphological searches. More recently, however, 
major progress is being made using the redshift information 
in the SDSS to pre-select systems with two discordant red- 
shifts. The major surveys to exploit this arc the OLS survey 
(Wilhs et al. 2006) and the SLAGS survey (Bolton et al. 
2006; Treu et al. 2006; Koopmans et al. 2006) which have 
already found 24 systems between them and axe Ukely to 
discover between 50 and 100 lens systems. These are nearly 
all radio-quiet systems with extended optical background 
sources which are useful for lens modelling. 

In this paper we address the question of making sur- 
veys for Icnscd radio sources more efficient where detailed 
redshift information is not available. There arc many sur- 
veys of objects, such as the FIRST survey, which in princi- 
ple contain several thousand objects in which the effects of 
lensing should be visible, but in which a direct search of all 
images at high enough resolution to detect lenses would be 
prohibitively expensive. This is the case even with the next 
generation of telescopes such as EVLA and e-MERLIN. To 
illustrate the extent of the problem, consider a survey of 
faint compact flat-spectrum radio objects. All flat-spectrum 
objects above 30mJy have already been searched for evi- 
dence of lensing in JVAS/CLASS, and to obtain ~10 times 
more objects, one would need to go a factor of 10 fainter. 
The increased sensitivity of future instruments allows for 
this factor of 10, but to get a factor of 10 more lenses, either 
10 times more targets need to bo observed (i.o a factor of 
ten more observing time is needed) or surveys will have to 
get a factor 10 more efficient. 

In this paper we explore ways of increasing efficiency of 
surveys for lens systems targeted at compact radio sources. 
The method attempts to make optimum use of existing 
general-purpose surveys. We describe the results of a small 
pilot of an "efficient" survey, calculate the expected lensing 
yield in a realistic situation, and consider the prospects for 
larger surveys of this type in the future. 



2 OBSERVATIONAL METHODS FOR 
FINDING RADIO LENSES 

2.1 Outline of the method 

As inputs we start with an optical and a radio survey of a 
patch of sky, neither of which have the resolution on their 
own to find gravitational lens system. For the purposes of 
the current investigation we distinguish five types of lens 
system: 

(i) Compact radio sources with images visible on the op- 
tical survey but with lens galaxies fainter than the optical 
survey limit 

(ii) Compact radio sources with images below the optical 
survey limit but with lens galaxies brighter than the optical 
survey limit 

(iii) Compact radio sources where both the optical images 
and the lensing galaxies are visible on the optical survey 

(iv) Moderately compact, but stcop-spoctrum, isolated 
radio sources with lens galaxies visible on the optical survey 

(v) Extended radio lobes associated with multiple 
sources, lensed by galaxies visible on the optical survey 



We note that the Haarsma et al. (2005) survey consti- 
tutes an attempt to find lenses of the fifth type, and there- 
fore do not consider them further. In cases where the ra- 
dio source is compact we also need to consider the possi- 
bility separately of detecting 2-imagc and 4-imagc systems. 
We also note that, if the SDSS is used as the optical sur- 
vey, JVAS/CLASS lenses would be classified as a mixture 
of types (ii) and (iii). Type (iii) cases are dominated by 
quasars from the strong-source end of the survey; at the 
fainter end, radio sources with flux densities <100 mjy are 
optically much fainter. In these cases, finding redshifts for 
them can be difficult even with large telescopes (Marlow et 
al. 2000; Falco, Kochanck & Mufioz 1998). 

There are two potentially useful efl'ects that may be 
seen, despite limited resolution in existing radio and op- 
tical catalogues, due to the fact that the radio catalogues 
are more sensitive to the source brightness than the lens 
brightness. First, lensing involves a displacement of the cen- 
troid of bright images away from the lens. This is particu- 
larly marked in the case of double-image systems, where the 
brighter image is of order 1.5-2 Einstein radii from the lens, 
but is also true for most four-image systems where the pair 
of merging images, which dominates the source brightness, 
is located one Einstein radius from the lens. We refer to 
this effect as the "optical-radio offset" (Fig. 1). In the most 
optimistic case, assuming that the optical image is domi- 
nated by light from the lens and not the lensed images, we 
can detect this offset at levels well below the resolution of 
both surveys. For instance, if FIRST is used as the radio 
survey, displacement can be detected at the 3-<t level down 
to separations of approximately 1" despite the 5" resolution 
of FIRST. A test based on the optical-radio offset becomes 
more powerful the more asymmetric the lens system. 

Second, the stretching introduced by the geometry of 
the multiple imaging can be detected even if the image 
separation is below the formal resolution, given sufficient 
signal-to-noise ratio. In the case of double-image systems, 
the stretching is parallel to the vector joining the centroid 
of the radio and optical source positions, and in the case 
of quad lenses it is perpendicular to it. We refer to this as 
the "misalignment angle" . Conveniently, at least for double- 
image lens systems, a test based on the misalignment angle 
becomes more powerful the less asymmetric the system. 

Figure 1 illustrates the overall geometry of these two 
situations. The misalignment angle M is defined as the an- 
gle between the radio component major axis and the line 
joining the radio and optical centroids, and should be close 
to 0° for double lenses and close to 90° for quads. Thus 
lensed systems should be predominantly located in certain 
defined regions of the M/optical-radio offset plane. This is 
the foundation of the method we are proposing. It should 
be emphasised, however, that although this method at first 
glance gives powerful detection of lens systems with a mini- 
mum of false positives, Figure 1 represents an ideal situation. 
When there are systematic astrometric errors within the sur- 
veys, or for type (iv) lenses with some extended structure, 
the situation is more complicated and it is important to use 
simulations and pilot observations to investigate the effec- 
tiveness of such methods in more detail. As a preliminary 
illustration that this method can work, we plot in Figure 
2 some actual CLASS data on three gravitational lens sys- 
tems. It can be seen that, provided the lensed images are not 



Efficient gravitational lens surveys 3 



FIRST 




Double image lens 

Figure 1. Diagram of the expected configuration for quad and 
double lens systems. In the case of double images, the FIRST 
centroid (dotted) is expected to be offset from the lens galaxy 
and to be extended roughly along the direction of the optical- 
radio separation due to the influence of the radio secondary (i.e. 
M=0). For quad lenses, the FIRST radio and SDSS optical galaxy 
will also be offset, because the close pair of images are very bright, 
and the radio axis will in general be nearly orthogonal to the 
optical-radio offset. 

optically dominant, the optical-radio offset can be used to 
distinguish cases of lonsing. In practice, the chances of the 
lensed object being optically bright diminish hugely with de- 
creasing radio flux density (Falco, Kochanek & Mufioz 1998; 
Marlow et al. 2000) and thus the test should become more 
efficient for fainter radio sources. The test is also particularly 
helpful for very asymmetric lenses, which are the best can- 
didates for searching for central images to give constraints 
on the mass distribution in the centre of the lensing galaxy 
(e.g. Boyce et al. 2006). 

2.2 The method applied to the FIRST and SDSS 
surveys 

Currently, the large-area surveys with the highest resolution 
are the optical SDSS (York et al. 2000) which covers an area 
of 7000 square degrees of the North Galactic Cap at an an- 
gular resolution of ~ l'.'4, and the FIRST radio survey at 
1.4 GHz (Becker, White & Helfand 1995) which covers a 
similar area at a resolution of 5". Both surveys claim astro- 
metric accuracy of a small fraction of an arcsecond. Neither 
survey has the resolution unambiguously to distinguish im- 
age splittings in average galaxy-lens systems of about 1". 
However, using a test based on the optical-radio offset and 
the misalignment angle should in this case be able to detect 
lens systems with separation of less than 1", given the astro- 
metric accuracy claimed by both surveys and the resolution 
of the FIRST survey of slightly less than 5". In principle, 
the test can be made more efficient for double-image lens 
systems by selecting misalignment angles close to 0° and for 
four- image systems by selecting angles close to 90° . 

A few JVAS/CLASS lenses would not be detectable on 
the SDSS. The SDSS flux limit is r'=23.1 (York et al. 2000) 



which corresponds to V ^ 22.6 (Jester et al. 2005), or to 
/ ~ 20.2 given a typical V — I colour of an elliptical galaxy 
at z ~ 0.5 (Bruzual 1983). Lens galaxy magnitudes are avail- 
able for all JVAS/CLASS lenses, and six are well below this 
limit: B0128+437 (Biggs et al. 2004), B0739-f366 (Marlow et 
al. 2001), B1127-I-385 (Koopmans et al. 1999), B1359-M54 
(Rusin et al. 2001) B1555-I-375 (Marlow et al. 1999) and 
B1938-I-666 (King et al. 1998). B1555-I-375 is in the SDSS 
region, and indeed is not detected. Therefore there is little 
prospect of using the methods outlined here for the detection 
of about 30% of radio source lenses. 

There are other useful radio surveys which can be used 
to establish spectral index of radio sources, which indicates 
whether the source is extended (steep-spectrum) or compact 
(flat-spectrum), in particular the Westerbork Northern Sky 
Survey (Rengelink et al. 1997) at 325MHz and the GB6 
survey (Condon et al. 1998) at 5GHz. 



3 A TEST OF THE SURVEY CONCEPT 
3.1 Sample selection 

Sources were selected from the regions of overlap between 
Data Release 3 of the SDSS and the FIRST survey. To avoid 
selection of radio lobes or more extended radio sources, for 
which the procedure outlined in section 2 probably does not 
work, only sources with fitted deconvolved sizes of 2'.' 5 or 
less from the FIRST survey were selected, and 3333 sources 
already observed by the CLASS survey were also rejected. 
An additional spectral index selection was achieved by re- 
quiring that all selected sources have FIRST integrated flux 
densities >9 mjy and have spectral index flatter than —0.5 
(S^ oc v") between the 1.4-GHz FIRST survey and detec- 
tions in the 325-MHz Westerbork Northern Sky Survey. The 
spectral index selection allows one to reject small compact 
steep spectrum sources with twin radio lobes lying at dif- 
ferent distances from the radio core which would otherwise 
contaminate the sample to am unacceptable extent. This pro- 
duces a sample of 3363 sources, of which 1831 have a poten- 
tial SDSS identification within 6". A further cut was per- 
formed based on the SDSS photometry in order to remove 
low-redshift z < 0.2 galaxies which are unlikely to be lenses 
by requiring that the fiuxes in the SDSS u, g and z bands 
obey the relation u — g < g — z. This in effect discriminates 
against galaxies which have a 4000-A break between the u 
and g bands. 1380 objects remain after applying this cut. 

In Figure 3 we plot the optical-radio offset versus mis- 
alignment angle M for the sample of 1380 candidates. In 
addition, we have simulated FIRST observations of the grav- 
itational lenses from CLASS, reducing their total fiux den- 
sity to 20 mJy and using their simulated FIRST character- 
istics, together with the known lens position, to calculate 
the optical-radio offset and misalignment angle M. Note 
that this assumes that the optical emission in the simulated 
lenses comes only from the lens galaxy, which is not true for 
the CLASS sample but is justified here because the lenses 
we are trying to find come from the <100 mJy population 
in which the sources are not optically bright. 

A final observing sample of 160 sources has been con- 
structed by choosing all objects with optical-radio offset be- 
tween 0'.'4 and 2", and with misalignment angles 0° < M < 
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Figure 2. Actual data for three CLASS lenses: B0850+054 (left), B1600+434 (centre) and B1030+074 (right), mapped using a restoring 
beam similar to the B-array 8.4-GHz observations used in this paper. The FIRST and SDSS positions are marked with a star and circle, 
respectively. Note that in the first two cases, the SDSS position is offset from the radio peak, reflecting the fact that the galaxy is offset 
from the brightest lensed image as in Figure 1. Note that this method does not work for B 1030+074 as here the lensed image is very 
bright in the optical. However, we argue that for significantly fainter radio sources, the likelihood of this happening is much less. 
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Optical— radio offset 

Figure 3. Misalignment angle M vs. optical— radio offset, in arc- 
sec, for sources in the FIRST-SDSS overlap region (specks). The 
CLASS lenses have been rcscalcd to 20mjy and mapped using 
artificial VLA snapshots and the same noise levels as FIRST. 
Double lenses (triangles) lie mostly in a region of the M:optical- 
radio offset plane well away from the bulk of the normal source 
population, as do the quad lenses (circled triangles). The excep- 
tions are the small-sepajration lenses which cannot be detected in 
this survey. The boxes enclose regions selected for observation, 
which contain a much larger fraction of gravitational lenses than 
the general population. 



30° and 60° < M < 90° . This cuts out the vast majority of 
unlensed sources, while kcepiug significant numbers of lenses 
(12 out of 22 of the simulated lens population remain in the 
final sample). Of these 12, we would expect ~9 to have lens 
galaxies dotoctablo with the SDSS. This implies that in ob- 
servations of this sample we should find a fraction 9/22 of 
lenses, while carrying out 160 observations instead of 3363, 
an increase in efficiency of nearly a factor 10. More lenses 
could be discovered at only slightly lower efficiency by in- 
cluding the 30° < M < 60° range of misalignment angle. 
This means that in principle one lens could be discovered 
per 60 sources instead of per 600, resulting in a lensing rate 
from an EVLA observation of 1 lens per 2 hours instead of 



the CLASS rate which was more like one lens per day from 
the VLA^. This factor of -^10 improvement is of the required 
order for the future extraction of 10 times larger samples of 
lenses. 



4 OBSERVATIONS AND RESULTS 
4.1 Observations 

4.1.1 VLA 

A total of 4.5 hours of VLA time was allocated to the pilot 
project in order to test the survey concept. Sources were ran- 
domly chosen from this sample according to available slots in 
the observing schedule. Three epochs of observations were 
obtained, 2005 March 03 (31 sources), 2005 March 17 (35 
sources) and 2005 March 29 (35 sources) for between 110 
and 120 seconds per source. The VLA in B-configuration 
was used at a frequency of 8.415 GHz and a total band- 
width of 100 MHz (2 IFs each of 50 MHz bandwidth). This 
arrangement gives a maximum baseline of 12 km, a nomi- 
nal resolution of ~ 0'.'7 and an rms noise level of approxi- 
mately 150 /iJybeam"^. Phase calibrators chosen from the 
JVAS survey were observed approximately every 20 minutes. 
These phase calibrators have a nominal astrometric accuracy 
of 40 mas or better. The phase stability throughout the ob- 
servations was good, typically with atmospheric phase drifts 
of < 30° over the whole observing session, and the accuracy 
of positions in the derived maps should therefore be better 
than about 100 mas. 

Unfortunately, the epoch from 2005 March 03 was af- 
fected by apparently random phase jumps between the 
Stokes R and L channels, mainly in the data on the target 

^ Wc note that if wc had considered sources selected by the 
method outlined here, but which had previously been observed 
with CLASS, we would have ended up with an observing list of 
69 sources and "discovered" one lens, CLASS B1600+434. The 
other lens in the appropriate region of sky, CLASS B1555-(-375, 
would have just been missed due to its lensing galaxy lying below 
the SDSS threshold 
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sources. Although correlated phases are obviously present, 
standard calibration and imaging procedures were unable 
to recover images of adequate quality for this epoch. The 
phases could be lined up by initial self-calibration to a point- 
source model, but at the risk of imposing spurious structure 
on the data and this epoch was therefore not used in the 
subsequent analysis. The total observed sample was there- 
fore only 70 sources, and the expected yield was, therefore, 
one gravitational lens system. 

The two good epochs were processed in a standard 
way with the NRAO Astronomical Image Processing System 
(aips). The flux density scale was normalised to an observa- 
tion of 3C286 in each case (Baars ct al. 1977) and an initial 
phase solution was constructed by interpolation across the 
phase calibrators. The amplitude and phase-calibrated data 
were mapped by hand with the difmap package (Shepherd 
1997) including CLEANing and one iteration of phase self- 
calibration for sources stronger than ~10 mjy. Final maps 
have been produced using uniform weighting and a restoring 
beam of 0f.'5, slightly below the nominal resolution of 0'.'67. 

4.1.2 MERLIN 

Two sources whose structures were ambiguous in the 
VLA observations, J1121-I-444 and J1316-I-575, were ob- 
served with the Multi-Element Radio-Linked Interferome- 
ter (MERLIN, Thomasson 1986) at an observing frequency 
of 4.994 GHz, with a 32-MHz bandwidth and a resolution 
of 50 mas. Observations were obtained on 2006 July 14- 
16 for 32h (J1316-I-575) and 17h (J1121-I-444) using four 
telescopes of the array (.JodrcU Mk2, Defford, Knockin and 
Cambridge). Phase calibrators were observed for 2.5 minutes 
in every 10 minutes, these being J1307+562 and J1115-I-416 
for the target sources J1316-I-575 and J1121-I-444. The point- 
source calibrator OQ208 and flux calibrator 3C286 were also 
observed, for 75 and 27 minutes respectively. Data were pro- 
cessed, and amplitude calibration applied, in the standard 
manner using the d-programme software available from the 
MERLIN site, and the flux density scale was set using the 
Baars et al. (1977) flux density for 3C286 of 7.3 Jy at 5 GHz, 
reduced by 30% to allow for resolution effects on the short- 
est available baseline (Knockin-Mk2). This reduction gives 
a flux density of 2.42 Jy for OQ208. Phase solutions were 
produced in aips for each scan of the two phase calibrators 
and interpolated to produce an overall atmospheric phase 
calibration. Maps were then made, using natural weighting 
in AIPS IMAGR and without self-calibration steps, yielding 
images with an rms noise of ~190 /^Jy/beam for J1121-I-444 
and ~110 /uJy/beam for J1316-I-575. This process produces 
clear detections of the phase calibrators. However, the maps 
of the two target sources do not show any detectable radio 
emission. 

4.2 Notes on resolved sources 

Maps of the resolved objects are presented in Figures 4 and 
5, and we discuss each case in turn. 

J1055-I-568. This source has a steep spectrum be- 
tween 1.4 and 8.4 GHz, declining from a peak bright- 
ness of 12mJy/beam in the FIRST survey to a peak of 
1.3mJy/beam in these maps. This is almost certainly be- 
cause at least one of the components is heavily resolved at 



8.4 GHz. Assuming the SDSS astrometry to be correct, a 
gravitational lens hypothesis is difficult to sustain as both 
radio components are on the same side of the optical ID, 
and the stronger component is close to the SDSS position. 
This source is almost certainly a steep-spectrum source with 
a radio jet. 

J1121-|-444. This source is still a candidate gravita- 
tional lens from the VLA observations. It is not a very good 
candidate, as the FIRST position is relatively far from the 
brighter 8.4-GHz component, suggesting a spectral index 
gradient across the source which is inconsistent with lensing. 
This impression is conflrmed by the MERLIN map (Fig. 5), 
which resolves out all radio emission. We conclude that both 
components are therefore extended on scales larger than the 
MERLIN resolution of 50 mas, and this is therefore not a 
gravitationally lensed compact source. 

J1220-I-476. The source is quite complex and appears 
to have a third component which is just detected in the radio 
map. Again, the components are almost certainly heavily 
resolved and the stronger radio component is close to the 
SDSS galaxy. It is thus difficult to claim this as a serious 
candidate for gravitational lensing. 

J1239-|-447. At flrst glance this object looks promis- 
ing, consisting of two apparently unresolved sources, one of 
which lies on the FIRST identification and one on the SDSS. 
The SDSS identification is with a. z = 2.04 object of r mag- 
nitude 19.8, and the overall radio spectrum is relatively flat 
(Sl.l = 0.48). A lensing identiflcation is however most un- 
likely, for the following reason. The fact that the FIRST 
centroid lies on one of the 8.4-GHz components, both of 
which have almost equal brightness, implies that the two 
components have different spectral indices. Such a differ- 
ence automatically rules out the two components as images 
in a gravitational lens system. Moreover, the SDSS quasar 
coincides with one 8.4 GHz component suggesting that the 
likely explanation for this system is that it is a one-sided jet 
in a radio quasar, where the core of the quasar is coincident 
with the SDSS identiflcation and the steep-spectrum radio 
jet is shifting the FIRST centroid of emission. 

J1259-I-519. As with J1239-I-447, this object is almost 
certainly a source with a one-sided jet. The FIRST identi- 
fication lies almost on top of one 8.4-GHz component, and 
the SDSS identification (a stellar object with r=21.5) on the 
other. In this system the overall spectral index is very steep, 
suggesting that at least one component is a heavily resolved 
radio lobe. 

J1316-I-575. This source is still a candidate from the 
VLA map, although the displacement of the 8.4-GHz compo- 
nent from the FIRST position and the position of the SDSS 
identification close to the brighter 8.4-GHz component make 
it an unlikely one. Again, however, the MERLIN map (Fig. 
5) resolves out both components and this is therefore not a 
gravitationally lensed compact source. 

J1321-I-476. The brighter component in the 8.4-GHz 
map is obviously resolved. In principle this might be thought 
to be the merging two images of a 4- image lens system. How- 
ever, the SDSS identification is in the wrong place either for 
a lensing galaxy or to be identified as a lensed radio image. 
Again, this system almost certainly consists of a compact 
core coincident with the SDSS position and a steep-spectrum 
radio jet which appears in FIRST. 
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distribution of FIRST controids reaches ~ 3a, or where the 
SDSS astrometry is in error. This is entirely consistent with 
the distribution in Figure 6, which shows that the distribu- 
tion of offsets between VLA and SDSS positions has a long 
tail, consisting of about one-third of the observed objects 
which are probably chance coincidences between FIRST and 
SDSS objects. 
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Figure 6. Position offsets for unresolved sources between our 
VLA measurements and the positions recorded in the FIRST and 
SDSS catalogues. 



J 1325+ 540. This radio source has a marginally de- 
tected secondary, but both components are very steep- 
spectrum and heavily resolved. 

J1441+567. This is potentially a good lens candidate. 
Unlike in the other cases, the FIRST position is between 
the two 8.4-GHz components detected on our radio map, 
suggesting that the two 8.4-GHz components have the same 
radio spectrum. However, the source has a relatively steep 
radio spectrum and there is a faint bridge of radio emission 
connecting the two components. This bridge is inconsistent 
with the results of 2-image lensing and entirely consistent it 
being a twin-jetted double radio source, whose radio emis- 
sion is not quite steep-spectrum enough to cause the source 
to be removed by the spectral index criterion. 

The conclusion is that no clear examples of gravitational 
lensing have been discovered in these observations. 



4.3 Nature of the observed sources 

In the last section we considered the sources in which our 
VLA observations revealed resolved structure. We now con- 
sider the sources that are unresolved by the 8.4-GHz VLA 
observations. Figure 6 shows the position offsets determined 
between our 8.4-GHz radio component and the FIRST and 
SDSS surveys. It can be seen that in the majority of the 
anomalous cases, the problem is an offset between the VLA 
position and the SDSS identification. 

The SDSS source density is approximately 0.0021 per 
square arcsecond. Hence in a search of 3764 sources within a 
radius of 2" we would expect 98±10 chance identifications. 
By reducing the sample of SDSS IDs from 2033 to 1538 us- 
ing the colour cut, we reduce the expected number of chance 
identifications to 74, and the selection on M and optical- 
radio offset reduces this further to 48±7 within the final 
sample of 160 (of which 70 were actually observed). There- 
fore about 30% of associations between FIRST sources and 
SDSS objects are likely to be chance coincidences, with the 
probability of a chance coincidence increasing as the SDSS 
object goes further from the radio source. The remaining 
70% are most probably cases in which there is an astrophys- 
ical reason for a optical-radio offset such as a core-jet struc- 
ture in a radio quasar, cases in which the tail of the error 



5 CONCLUSION 

We have shown that it is in principle possible to increase 
the efficiency of lens searches in compact radio sources by 
up to a factor of --^10, resulting in a lensing rate of 1 in 
60 rather than 1 in 600 in the most favourable case. This 
can be achieved by using sub-resolution position shifts and 
slight resolution effects within the radio FIRST survey and 
the SDSS survey, supplemented by radio observations. We 
present a preliminary programme with the VLA and MER- 
LIN in which we would have expected to detect 1 lens out of 
an observed sample of 70; we have eliminated all candidates 
and therefore detect no lenses. 

The major difficulty with this method is that its effi- 
ciency is fundamentally limited by the high density of optical 
objects at the SDSS r' = 23.1, which means that each true 
lens system will be accompanied by several tens of false- 
positives from random coincidences with unrelated SDSS 
objects. In principle, this could be reduced further by more 
stringent colour cuts to eliminate SDSS identifications which 
are unlikely lensing galaxies. Further false-positives are pro- 
vided by the tail of sources with 3-cr errors in the optical 
astrometry, and efficiency would be increased dramatically 
by a 30% reduction in astrometric errors. The final prob- 
lem is similar to the limiting factor in the CLASS survey, 
namely that one runs out of radio sources even at a lower flux 
level. In the current survey we have observed a third of the 
available radio sources with Si.4ghz > 9mJy in the SDSS 
footprint. Further progress would require a larger volume of 
sky observed by SDSS, or a deeper radio survey at a fre- 
quency sufficiently different from 1.4 GHz to allow spectral 
selection. Alternatively, one could relax the spectral index 
criterion in the sample selection, at a cost to efficiency. 

A major advantage with the method, however, is that 
lens systems are radio loud and can be confirmed by purely 
radio observations rather than requiring high-resolution op- 
tical observations using scarce resources such as the Hub- 
ble Space Telescope. Moreover, in the future this method 
may find more widespread applicability as all-sky surveys 
with astrometric precisions of 0.1 arcseconds or better be- 
come commonplace. Such surveys are likely to be done in the 
next ten years with, for example, the Large Synoptic Sur- 
vey Telescope (LSST, Claver et al. 2004) in the optical, and 
in the radio with the Low-Frequency Array Radiotelescope 
(LOFAR, Rottgcring ct al. 2005). These surveys should give 
rise to samples of many millions of radio sources with opti- 
cal identifications and spectral index information, allowing 
the efficient discovery of thousands of lenses with modest 
foUowup time. 
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Figure 4. Resolved sources from the VLA observations. The FIRST and SDSS positions are indicated by five-pointed stars and circles 
respectively. Top row, left to right: J1055+568, J1121-I-444, J1220+476. Middle row, left to right: J1239+447, J1259+519, J1316+575. 
Bottom row, left to right: J1321+476, J1325+540, J1441+567. 
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